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Abstract 

We find power-law behaviors of grain velocity in both propagation and 
backscattering in a gravitationally compacted granular chain with nonlin¬ 
ear contact force. We focus on the leading peak of the velocity signal which 
decreases in a power-law d~°‘, where d is the location of the peak, as the 
signal goes down. The ratio of backscattered to incident leading velocity also 
follows a power-law df^, where di is the depth of impurity. The up-going 
backscattered signal is nearly solitary. Therefore, the overall change of the 
leading velocity peak is given by the power-law d- We find a = 0.250 

and f3 = 0.167 for the Hertzian contact force. 
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Recently, physics of granular materials attracts great interest, since the materials are 
ubiquitous around us and their properties are unique and also useful in many applications 
1]^. The propagation of sound or elastic wave in granular medium is also one of interesting 
subjects related to the properties of granular matter [2], A rather simple system, the granular 
chain with nonlinear contact force including Hertzian contact |0], has been revived by hnding 
a soliton in transmitting elastic impulse. This soliton found in a loaded horizontal Hertzian 
chain was hrst predicted by Nesterenko @] and its experimental verihcation was performed by 
Lazaridi and Nesterenko and recently by Coste et al. P] . Even though three-dimensional 
granular systems may not follow simple Hertzian contact force law due to geometrical effect 
0, the simple nonlinear contact law is still interesting because the existence of solitary wave 
may provide a possibility to get information inside granular matter and the one-dimensional 
system usually provides the starting point for studying more complex systems. 

In addition to studying horizontal chain, the gravitationally compacted chain and column 
of grains have been studied by Sinkovits and Sen P-jl0[| in terms of molecular dynamics 
simulations. The two cases are not the same in the following sense that the former shows 
robust solitary wave within certain amount of loading, while the latter shows solitonlike 
wave whose waveform disperses little by little as it propagates down. Sen at el. 0 studied 
backscattering of the solitonlike wave by impurity in a gravitationally compacted chain to get 
the information on buried impurity such metal-poor land mines. They obtained interesting 
relations between the phase of grain-velocity signal and the mass of impurity. The result 
reflects the well-known phase-density relation in the connected string with different densities. 

In this Letter, we focus on the propagating and backscattering characteristics of the 
leading part of grain-velocity signal in the gravitationally compacted granular chain with 
Hertzian contact force especially. We also treat other types of nonlinear forces. We report 
interesting power-law behaviors of the velocity signal. This work is attractive because it 
supplies new laws of signal propagation in the gravitationally compacted granular chain 
with nonlinear contact force and a possibility of identifying the location of buried impurity 
by observing the peak-velocity of the leading part of returned signal at surface. 
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To study the dynamics of grains in the gravitationally compacted Hertzian chain, we 


solve numerically the equation of motion of a grain at Zi, 


mzi = -a[{Ao - (zi - - {Aq - (zi+i - Zi)}^^'^] 

- mg, 


( 1 ) 


where m is the mass of grain, Aq = Ri + Ri+i, and a is the constant dehned in Eq. (2) 
below. We neglect plastic deformation in Eq. (1). This equation of motion comes from the 
Hertzian interaction energy between neighboring granular spheres which is given by H 




•^ 2 , 2+1 


2 , 2+1 5 


( 2 ) 


5T^ \Ri Ri+i 

where is the overlap between two adjacent spherical grains and is given by = 

Ri + in which Ri is the radius of the sphere whose center is at position i and 

is the distance between i and z + 1, and 


D = 


3/1-+ 1--2 


+ 


a. 


i+l 


( 3 ) 


4 \ Ei Ej+i J 

where cxj, cTj+i and Ei, Ej+i are Poisson’s ratios and Young’s moduli of the bodies at neigh¬ 
boring positions, respectively 0]. The equilibrium condition 


N 


bci 3/2 

9 2^ mj = — 5^+1 
j=i+i ^ 


( 4 ) 


has been used for the z-th grain of a chain of N grains. A criterion for the initial impulse 
exists to make Eqs. (1) and (2) valid P,pdj|. Therefore, in this work, we choose comparatively 
weak impulses to satisfy this criterion. 

Our system is a vertical chain of A = 10^ grains which has hve consecutive impurity 
grains which are different only in mass from other spherical grains. The mass of impurity 
is taken as half of the mass of medium grains. As a calculational tool, we use the third- 
order Gear predictor-corrector algorithm []^ and the same program units as those used in 
Ref. [lOT. That is, the units of distance, mass, and time are 10“^m, 2.36 x 10“®kg, and 
1.0102 X 10 ^s, respectively. These units gives the gravitational acceleration g = 1. We set 
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the grain diameter 100, mass 1, and the constant a of Eq. (1) 5657 for molecular dynamics 
simulation. 

We focus on the peak value of the very leading velocity signals from the moment of 
initial impulse imposed on the hrst grain at top of chain to the moment of return to the 
hrst grain. Figure 1(a) shows the snap shots of propagating velocity signal passing different 
depths. Initial impulse velocity 0.1 in program units has been applied to the hrst grain. 
One can easily see from Fig. 1(a) that more grains are involved in the signal, i.e., signal 
disperses as time passes. We hnd that even though loading on the grain increases in depth 
due to gravity, total kinetic energy and potential energy of the grains involved in the signal 
remains constant independently. Therefore, the decrease of peak velocity results from only 
dispersion of the signal. Fig. 1(a) shows that the propagating signal is not solitary in the 
gravitationally compacted chain. This is quite different from the case of loaded horizontal 
chain [^,|^. 

In this work, we are interested in the effect of gravity which changes loading of chain 
continuously and the time- or depth-dependent behavior of the leading peak shown in Fig. 
1(a). Fig. 1(b) is the logio-logio plot of the leading velocity peak versus its location. A 
clear power-law behavior is seen in Fig. 1(b). The explicit expression for the leading peak 
velocity is given by 


v{d) 

Vi 


Adr 


(5) 


where Vi is the initial velocity, d is the distance from top to the leading peak, and A is the 
intersection which is not equal to one because of boundary effect at the top of chain. But 
the signal quickly goes to the stable form shown in Fig. 1(a). We will show this in Fig. 2 
. Figure 1(b) obtained for Hertzian contact gives rise to the power-law index a = 0.250. 
Figure 2 drawn in logio-logio plot shows the propagating behaviors of the leading peak of 
the velocity signal for various nonlinear contact forces, such as n = 5/2, 3, 7/2, and 4 of 
the potential V{6i^i+i) = The abscissa of Fig. 2 denotes elapsed time after impulse. 

One can see that the initial impulse quickly follows the power-law of Eq. (^ in 0.1 second 


4 



which corresponds to passing 40th grain in length. Interesting things in Fig. 2 are (1) the 
slopes are very close to each other even though there is crossing, therefore, one expects 
similar propagating behaviors for these nonlinear contact forces, and (2) Hertzian case is 
well-separated from other nonlinear cases. It is unclear right now why Hertzian is special. 

The signal performs backscattering when it meets impurity. Figure 3(a) shows typical 
backscattering by five consecutive impurity grains of mass 0.5 and depth 500 in program 
units. The total energy of the incident signal is divided into transmitted and reflected parts, 
i.e., ~ + '^hack- The backscattered leading peak diminishes abruptly compared 

with second leading peak in the process of scattering. We call the ratio of backscattered 
to incident leading peak the reflection coefficient of the leading peak velocity for present 
analysis. We find a remarkable phenomenon that the reflection coefficient of the leading 
peak velocity also follows a power-law in depth of impurity. The reflection coefficient, of 
course, depends on the mass of impurity. We restrict this study to the case of impurity mass 
0.5 while the mass of medium grain is 1. Fig. 3(b) is the logio-logio plot of the reflection 
coefficient, i.e., the ratio of backscattered leading peak to incident one versus the location 
of impurity dj. A clear power-law is seen in Fig. 3(b) and the expression is given by 




( 6 ) 


where v^^di) is the leading peak velocity after backscattering by impurities at depth d* and 
B is the intersection of logio-logio plot. Figure 3(b) obtained for Hertzian contact gives rise 
to the power-law index f3 = 0.167 for the impurity mass 0.5. We take the value of Vf,(dj) 
at the moment of separation from the tail of incident signal. We confirm that conservation 
laws of energy and momentum are always retained. 

This power-law behavior is an unexpected result. We thought this may stem from the 
change of loading as the location of impurity changes in the gravitationally compacted chain. 
We checked this idea for the horizontal chain by changing the amount of loading. We found 
that there is no change in the reflection coefficient upon loading for the solitary wave in the 
horizontal chain. We also confirm that the amplitude of incident solitary wave in horizontal 
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chain cannot change the reflection coefficient either. Therefore, we conclude that the power- 
law behavior in Eq. is a characteristic of nonsolitariness of the wave in gravitationally 
compacted granular chain. 

We also check the propagation of backscattered wave. It is interesting to see that the 
backscattered waveform is very solid and it shows negligible change in the height of the 
leading peak within the depth we considered in this work. Also interesting thing is that 
the propagating behavior of the backscattered velocity signal is not power-law but linear in 
depth as shown in Fig. 4 drawn in absolute scale. The saw tooth form stems from finite 
radius of grain. The trend of changing leading peak in Fig. 4 is flat in the velocity scale 
used in this analysis. We And that the number of grains participating in the signal increases 
in the case of signal going down, while it remains nearly constant in the up-going situation. 
Therefore, there is little dispersion when the signal goes up in the gravitationally compacted 
chain. In addition, no kinetic and potential energy change either in this case. The linear- 
law of Fig. 4 does not play a role for the shallow region but play an appreciable role when 
backscattering produces at a deep enough place. 

As a summary of the above analysis, we obtain the overall power-law relation between 
returned velocity Vr at top and the depth of impurity dj. Since the kinetic and contact 
potential energy is transferred to the kinetic and gravitational potential energy of the first 
grain at top, the leading backscattered peak velocity at grain 2, i.e., Vh{2) ^ —0.004 in Fig. 
3(a), is quite different from the leading peak velocity at grain 1 which we define the returned 
leading peak velocity Vr- Figure 5 shows the logio-logio plots of Vr/vi versus the location 
of impurity di for different contact forces and initial impulse velocities Vi = 0.1,0.05, and 
0.01 in program units. One can see a big difference between Vr and nb(2). We separate 
artificially the lines of Fig. 5 into two groups by dividing the data of V{6i^i+i ) = <,+. by 
1.1 to avoid overlap with those of Hertzian interaction One can see that 

the power-law is independent of initial impulse. The equation of the straight lines of Fig. 5 
is written as 
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V, 


(7) 


Vi 




where the coefficient R{vi), the intersection of the logio-logio plot, contains the energy 
transfer at top and the coefficients A and B of Eqs. (5) and (6), respectively. R{vi)^s weakly 
depends on the amonnt of initial impnlse and the type of contact force law as shown in Fig. 
5. The power-law index a + f3 = 0.417 for Hertzian contact and 0.378 for another. From the 
analysis above, we know that the index is composed of two different sonrces, i.e., the index 
for downward propagation a and the index for backscattering [3. One can conhrm this by 
adding the slopes of Fig. 1(b) and Fig. 2(b) for Hertzian chain. 

In conclusion, we solve numerically the many-body equations of motion of the gravi¬ 
tationally compacted granular chain coupled by nonlinear contact forces and analyze the 
propagating behaviors of the leading peak of the velocity signal and the backscattering 
characteristic by buried impurity. We hnd two power-law relations in Eqs. (|^) and (|]) for 
propagation and backscattering, respectively. The former contains only the effect of dis¬ 
persion of nonsolitary wave as it propagates. The latter, however, has nothing to do with 
dispersion but the phenomenon occurring at the point of scattering. Comparing with the 
case of solitary wave in loaded horizontal chain, we conclude that Eq. (j^) is a special feature 
of nonsolitary property of the gravitationally compacted granular chain. The backscattered 
up-going wave, on the other hand, is a nearly solitary wave and does not follow power-law but 
follow very weak linear-law unlike the down-going wave. We conclude that this interesting 
feature is the characteristic of the gravitationally compacted granular chain with nonlinear 
contact force. The solitary behavior of the up-going wave in the gravitationally compacted 
granular chain may be useful to hnd buried impurity like nonmetallic land mine. 

InEq. (0), we obtain the combined power-law on the location of impurity for the returned 
leading peak velocity at top of chain. This power-law also exists for other nonlinear contact 
forces as well as Hertzian as we show in Fig. 5. We therefore conclude that above features 
and power-laws are generic for any nonlinear contact forces of type F{5i^i+i) oc ^^le 

gravitationally compacted granular chain. 
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Figure Captions 


Fig. 1 ; (a) Snap shots of propagating wave forms in a gravitationally compacted granular 
chain with Hertzian contact law. (b) Logio-logio plot of v{d) vs. depth d for the 
leading peak in (a). Slope of the straight line is —0.250. 

Fig. 2 : Logio-logio plot of n((i) vs. d for the leading peak for various interaction potentials, 
i.e., n = 5/2, 3, 7/2, and 4 of H(dj,j+i) = od/j+i- 

Fig. 3 ; (a) Snap shot of scattering by hve light impurities at 500. (b) Logio-logio plot of 
Vb{di)/v{di) vs. di- Slope of the straight line is —0.167. 

Fig. 4 ; Plot of leading peak vs. its location for the backscattered part of the velocity 
signal. Absolute scales are used for both abscissa and ordinate. Five impurities are at 
700. 

Fig. 5 : Logio-logio plot of Vr/vi vs. di for the leading peak for various initial velocities, 
Vi = 0.1,0.05,0.01. Upper three lines are Hertzian (n = 5/2) cases. Lower three lines 
(n = 3) are lowered artihcially. Slopes are 0.416 (n = 5/2) and 0.378 (n = 3). 
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